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Results are presented on a partial wave analysis of the lot] final state produced in n~p interactions 
at 18 GeV/c where uj — > ir + ir~ir°, ir° — > 2j , and r\ — > 27. We observe the previously unreported 
decay mode oj(1650) — > lot) and a new l" 1 meson state /ii(1595) with a mass M = 1594(±15)(+g°) 
MeV/c 2 and a width V = 384(±60) ( + 1 7 ° ) MeV/c 2 . The /n(1595) state exhibits resonant-like phase 
motion relative to the tj(1650). 



INTRODUCTION 



Studies of meson spectra via strong decays of 
hadrons provide insight regarding QCD at the confine- 
ment scale. These studies have led to phenomenological 
models such as the constituent quark model. However, 
QCD demands a much richer spectrum of meson states 
which includes extra states such as hybrids (qqg), mul- 
tiquarks (qqqq), and glueballs (gg or ggg). Experiment 
E852 at Brookhaven National Laboratory is an experi- 
ment in meson spectroscopy configured to detect both 
neutral and charged final meson states of ir~p collisions 
in a search for meson states beyond those compatible 
with the constituent quark model. 

The apparatus was located at the Multi-Particle Spec- 
trometer (MPS) of Brookhaven's Alternating Gradient 
Synchrotron (AGS). The AGS delivered an 18 GeV/c ir~ 
beam to a fixed liquid hydrogen target at the MPS. The 
MPS facility was augmented with additional detectors 
designed specifically for E852 which consisted of 3 inte- 
gral regions: target, charged tracking, and downstream 
regions(see Figure Q). 

The target region was located in the middle of the MPS 
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FIG. 1. The plan view of E852 Apparatus located at 
the Multi-Particle Spectrometer of Brookhaven's Alternating 
Gradient Synchrotron. 



dipole magnet (1 Tesla) and contained the following el- 
ements: a 30.5 cm long liquid hydrogen target, a four- 
layer cylindrical drift-proportional wire chamber jl] , and 
a 198 block CsI(Ti) barrel veto detector @. 

The downstream half of the MPS magnet housed the 
main components of the charged tracking region. It con- 
sisted of 3 proportional wire chambers (TPX1-3) and 6 
drift chamber modules(DMl-6) each with seven-layers. 
In addition to the tracking chambers, there were two scin- 
tillation veto counters (CPVB and CPVC) and a window- 
frame lead-scintillator sandwich veto detector (DEA). 

The downstream region contained a 3000 element lead- 
glass calorimeter (LGD) [|| for detecting and measuring 
the energy of the forward going gammas and a large drift 
chamber (TDX4) located directly in front of the LGD 
for tagging charged particles entering the LGD. More 
detailed descriptions and general features of the E852 
apparatus, data acquisition, event reconstruction and se- 
lection are given in References Q-^] . 

The tor] final state is of considerable interest because it 
has been virtually unexplored, and since it has not been 
observed in the decay modes of known mesons. Also, 
Close and Page 0, using an extension of the flux- tube 
model of Isgur and Paton || , suggest that the isoscalar 
qqg hybrid should decay dominantly to pn 

and uir\. 

The exclusive uirj system has been studied by two 
experiments previously. The GAMS Collaboration ob- 
served less than 100 tor] events produced in n~p inter- 
actions Q and claimed to observe a narrow (less than 
50 MeV/c 2 wide) structure at 1650 MeV/c 2 . Photopro- 
duction of nearly 100 uirj final state events was observed 
by the Omega Photon Collaboration |l(| . They reported 
observing a peak in the tor) spectrum at a mass of 1610 
MeV/c 2 but with a width of 230 MeV/c 2 . Both experi- 
ments suffered from having too little data and could not 
perform a partial wave analysis. We report here the re- 
sults of a partial wave analysis of approximately 20000 
exclusive ior\ events. 



FEATURES OF THE DATA 

During the 1995 AGS data run, a sample of 750 million 
triggers was acquired of which 108 million were of a type 
designed to enrich the yield from the reaction n~p — > 
ri7r + 7T - 47. About 6 million events containing 7r + 7r~47 
were fully reconstructed. 

The data were kinematically fitted to select events con- 
sistent with a n7r + 7r~7r°?7 hypothesis. To eliminate bro- 
ken and incorrectly hypothesized events, kinematic fits 
with confidence levels less than 5% were excluded. 

There are six ways to combine the 4 7's into pairs, 
and one problem in interpreting the results of the kine- 
matic fit is deciding between these 77 combinations for 
ambiguous solutions. For example, twenty-nine percent 



of the fitted rwr 4 ~ir~ tt c 'rj events were also found to fit the 
7i7r + 7r~7r 7r° hypothesis. Events of this type with a con- 
fidence level of at least 0.01% for the n7r + 7r~7r°7r° hy- 
pothesis were removed from the data sample. A total of 
113000 7i7r + 7r~7r°r/ events remained for further analysis. 

The two prominent resonances seen in the 7r + 7r~7r° in- 
variant mass distribution (see Figure ^|) are consistent 
with the well known 77 and ui. A Gaussian plus a second 
order polynomial fit to these peaks results in a measured 
mass and width of the 7/ of (550.4 ± 0.2MeV/c 2 ) and 
23.4 ± 0.2MeV/c 2 and of the u of (787.5 ± 0.3MeV/c 2 ) 
and 38.3 ± 0.3MeV/ c 2 respectively. The mass values are 
in good agreement with the PDG Q mass values for 
the 77 and lj, whereas the fitted widths are a measure of 
the experimental mass resolution of the E852 apparatus. 
The 77 — > 7T + 7r~7r° events have been selected for a study 
of the 7777 system which is currently in progress. 
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FIG. 2. The ir + ir~ir° effective mass distribution. Two 
prominent resonances are seen: the r\ and the u). Events with 
an invariant 7r + 7r _ 7r° mass in the region of 750MeV/c 2 < 
mass(37r) < 830A/eV/c 2 were used to select lu events. 

Events with an invariant ir + ir~ir° mass in region of 
750Mey/c 2 < mass(37r) < 830AfeV7c 2 were used to se- 
lect to events. Figure || exhibits the J p — l - nature of 
the to. Displayed is the to decay matrix element squared 
(A.) (12|: 



where 7? + and tt~ are the three-momentum vectors of 
the the ir + and n~ in the 3-7T rest-frame. The 3-7T decay 
amplitudes are constructed from the n momentum vec- 
tors, and due to the overall negative intrinsic parity of 
37r, a J p = 1~ decay amplitude has to be built out of a 



pseudovector combination q |T^] : 



q 



7f°xn 



(1) 



The A distribution should rise linearly for omega events, 
whereas it would be constant for 37r events distributed 
according to phase-space. The 37r background under the 
to is considerable, but the decay information of the uj can 
be used to weight oj events more highly than non-omega 
events (see the section on Background Study). 
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FIG. 3. Distribution of events as a function of the omega 
decay matrix element squared, A. Pure J p = 1~ events would 
exhibit a distribution which linearly increases with A whereas 
phase-space events should be flat. 

In Figure [|, the lot] invariant mass distribution of 19530 
lot] final state events is shown using 40 MeV/c 2 mass 
bins. The invariant mass distribution rises rapidly near 
threshold, then increases at a slower rate to a maximum 
near 1600 MeV/c 2 . But, in general the mass distribu- 
tion shows no clear evidence for structures indicating the 
presence of resonant states. 



PARTIAL WAVE ANALYSIS 

A partial wave analysis (PWA) of the data was carried 
out using the BNL PWA program. A general description 
of the BNL PWA program is given in Reference |h|] . The 
methods used in the PWA are based on the formalism 
of the Isobar Model jlij. For related material on this 
formalism sec References 

The data used in the PWA are shown in Figure ||. 
A PWA, taking into account both nucleon spin-flip and 
spin-nonflip contributions, was performed in fifteen 50 
MeV/c 2 mass bins from 1320 to 2070 MeV/c 2 . 
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FIG. 4. The uir/ effective mass distribution from the reac- 
tion ir~p — > Lorjn (not corrected for acceptance). 

An extensive series of tests were performed to judge 
the stability of the PWA fit. These include a systematic 
study of the PWA by varying the: 

• allowed waves in the fit, 

• fit starting parameters, 

• data selection cuts, and 

• mass bin size. 

In each test, the general features of the PWA fit remained 
unchanged. 

An appropriate minimal set of partial waves needed to 
describe the data was determined by performing PWA 
fits with varying sets of partial waves. Initially, all par- 
tial waves with L < 4 were introduced (L is the relative 
orbital angular momentum between the lo and the rf). 
Partial waves were discarded from the fit if their contri- 
butions were negligible and if the removal had little effect 
on the remaining features and goodness of the fit. It was 
determined that a minimum set of 11 waves, which in- 
cluded a non-interfering isotropic background wave, were 
required in order to achieve a reasonable agreement of the 
experimental data and the PWA fit results (see section 
on PWA fit quality). A notable feature was that the ex- 
otic partial waves for J PC = - ~ and 2 H were found 
not to be required in the analysis. 

Table [| lists the partial waves included in the final 
PWA fit. Also included was an isotropic, non-omega, 
non-interfering background wave. The amplitudes are 
expressed in the reflectivity basis which takes into ac- 
count parity conservation in the production process by 
a transformation of helicity states to eigenstates of the 



TABLE I. The minimal set of partial waves required in the 
PWA of the ojr) system. 



Partial Waves J PC M e L 
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production plane reflection operator fl7| . For irp interac- 
tions, the reflectivity coincides with the naturality of the 
exchange particle and amplitudes of different reflectivity 
e = ± do not interfere. 



QUALITY OF THE PWA FITS 

The best test of the PWA fit is to compare the events 
generated according to the PWA to the data itself. For 
this step, Monte Carlo(MC) events were subjected to a 
simulation of the E852 apparatus and weighted with the 
probability obtained from the PWA fit. These weighted 
MC events should mimic the data for successful PWA 
fits. Since the fits were performed in independent mass 
bins, the MC events were normalized to the number of ex- 
perimentally observed events in the corresponding mass 
bin: 
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FIG. 5. The quality of the fit is shown by a comparison 
of angular distributions for the data (with error bars) and 
PWA fit predicted Monte Carlo (dashed histogram) : the Got- 
tfried-Jackson frame a) Cos(#), and b) (p. 
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Figures |^ and |6| shows a comparison of the weighted 
MC events to the data for the: Gottfried- Jackson angles 
il, and helicity angles £lh- The weighted MC events is 
in good agreement with the data. A more stringent test 
compares the angular moments, H (I, m; L, M )f| . The 
results for this comparison (see Reference Q ) show good 
agreement between data and predicted events. Overall, 
the comparisons of the distributions demonstrate that 
the data is very well described by the PWA fit results. 
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FIG. 6. The quality of the fit is shown by a comparison of 
angular distributions for the data (with error bars) and PWA 
fit predicted Monte Carlo (dashed histogram): the omega he- 
licity frame c) Cos(9h), and d) 4>h- 



PWA RESULTS 

Figures and || show the acceptance corrected partial 
wave intensities. The most significant partial waves in 
the analysis are the l" 1 0+5* and 1~~1+P natural parity 
waves. These two waves also exhibit relative phase mo- 
tion indicative of resonant states (see Figure ^c). Note 
that since the amplitudes which interfere have either 
purely real or purely imaginary decay factors there is 
an unavoidable but trivial ambiguity in the sign of the 
relative phase. The remaining natural parity waves are 
small and do not exhibit clear phase motion. The con- 
tributions of the two unnatural parity exchange waves 
tend to become more important at masses higher than 
1700 MeV/c 2 . Their relative phase is not well deter- 
mined and exhibits large errors; therefore it cannot be 
reliably used to study resonance behavior. Nonetheless, 
the peak in the 2~~1~P intensity was fitted with a Brcit- 
Wigner resonance shape giving a mass of 1830±8 MeV/ c 2 
and a width of 86 ± 31MeV/c 2 . This is interesting since 
the quark model predicts J PC = 2~~ states, yet none 
have been observed. A prediction by Godfrey and Is- 
gur suggests a mass of 1700 McV/c 2 for the lowest lying 
J p c =2— qq state @. 
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FIG. 7. Acceptance-corrected partial waves intensities for 



(J^M'L): 
and l+'l+S. 



1-P,3—0~F, I"" 1+P, l+'O+S, 1+-Q+D, 



The only previously reported isoscalar J PC = 1 + 
states are the /ii(1170) and /ii(1380) Q. Both of these 



states have masses much lower than the observed struc- 
ture in the 1+-0+S wave. Also, the /ii(1380), believed 
to be the ideal mixed ss state, is not expected to be pro- 
duced in n~p interactions due to OZI suppression. For 
the isoscalar J PC = there are two states listed by 

the Particle Data Group (PDG) in the 1000-2000 MeV/c 2 
range: the w(1650) and the w(1420) Q. The w(1650), 
which has a PDG mass of 1649 ± 24 MeV/ c 2 and a width 
of 220 ±35 MeV/c 2 , is in good agreement with the mass 
and width of the observed structure in the 1" 1+P wave. 

In order to determine if the 1 H + S and the 1~~1+P 
observed structures are consistent with resonance behav- 
ior, a mass dependent analysis(MDA) of the PWA re- 
sults was performed. The input quantities included the 
l" 1 0+5* and 1~~1+P partial wave intensities and their 
relative phase for each mass bin. The errors were cal- 
culated using the error matrix from the PWA fit, which 
takes into account correlations between the intensities 
and relative phases. Relativistic Breit-Wigner forms 
were used to parameterize the amplitudes for the two 
waves. The parameters of the MDA fit included the 
Breit-Wigner masses, widths, and intensities. An ad- 
ditional constant parameter was included to allow for a 
relative constant production phase between the waves. A 
series of mass dependent fits were performed for different 
hypotheses and for the different phase solutions. 



Events: 1 ID 



Events: 2 " O P 



*5l500 



^1000 

g 500 

> 





>1500 
O1000 

1 500 

> 

m 

o 



1500 
1000 
500 




I 



1.4 1.6 1.8 2.0 GeV/c 2 1-4 1-6 1.8 2.0 

Events: 2 " 1 + P Events: 3" 1*F 



1500 
1000 
500 



u_ 



1.4 1.6 1.8 2.0 GeV/c 2 1-4 



1.6 



1.8 2.0 



Events: background 



*u 1500 

s 

g 1000 

g 500 

> 

w 

o l 
FIG. 8. 



1500 
1000 
500 




1.4 1.6 1.8 2.0 GeV/c 2 1.4 1.6 1.8 2.0 

Mass[cor| ] Mass[C0T| ] 

Acceptance-corrected partial waves intensities for 



(J PC M e L): 1 + -1 + D,2~0+P,2— 1 + P,3 __ 1+P, all 
and background wave. 



The primary goal of this analysis was to understand 
the structure in the l" 1 wave. A fit of the 1 in- 
tensity to a single Breit-Wigner resulted in a mass of 
1700(±20) MeV/c 2 and a width of 250(±50) MeV/c 2 for 
a x 2 / (Tdof) of 0.77. These values are in reasonable agree- 
ment with the PDG values of the w(1650), and since the 
w(1650) is rather well known (see references 0j), we used 
the 1 signal as a "calibration" by identifying it with 
the w(1650) and fixing the 1" Breit-Wigner parameters 
to the PDG values and then fitted the for the l" 1 pa- 
rameters. This identification would constitute the first 
observation of the ujij decay mode of the k>(1650). Un- 
der these assumptions, a simultaneous MDA fit to the 
F intensity and the F /1 — - relative phase resulted 
in a x 2 /(!6do/) of 0.71. The l+^O+S mass and width 
resulting from this fit were 1594(±15)(+g°) MeV/c 2 and 
384(±60)(+ 1 7 ° ) MeV/c 2 respectively. The quoted errors 
correspond to statistical and systematic uncertainties, re- 
spectively. The systematic errors were estimated by fit- 
ting the PWA results obtained for different set of partial 
waves which varied from 9 to 22 waves. Figure ^ displays 
the results of the mass dependent analysis for this fit: ^a, 
and show results of this fit overlayed on the partial 
wave intensities and relative phase. Figure fyl shows the 
absolute Breit-Wigner phases for 1" and F and the 
relative constant production phase. On the other hand, 
fixing the 1 parameters to the Breit-Wigner fitted val- 
ues (mass= 1700 MeV/c 2 and width = 250 MeV/c 2 ) and 
performing a similar MDA fit resulted in a x 2 / (lQdof) 
of 1.28 and a 1+^0+5 mass = 1600(±20) MeV/c 2 and 
T = 455(±80) MeV/c 2 . 

If the F + S object is interpreted as being caused by 
a single resonance, then this state, an /ii(1595), does not 
coincide with any known states. Possible interpretations 
for an /ii(1595) include an hi radial excitation, an hi 
hybrid, or a radial-hybrid mixture. In the Godfrey-Isgur 



potential model 18 1, a mass of 1780 MeV/c 2 is predicted 



for the 2 3 Pq hi radial excitation. This is approximately 
180 MeV/ c 2 higher than the value from our analysis. 
However this model also predicts a high mass for the 
/ii(1160) and a high mass for the /i x (1380) (1220 MeV/c 2 
and 1470 MeV/c 2 , respectively), therefore one might ex- 
pect the radial state to actually lie 100-200 MeV/c 2 below 
the Godfrey-Isgur prediction. 

A calculation by T. Barnes p0| (also see Reference 21 ) 



using a 3 Pq model suggests that a 1700 MeV/ c 2 hi ra- 
dial excitation should decay almost equally to F S and 
F D partial wavesfj] which is not consistent with what 
we observe. Alternatively, since the flux-tube model pre- 
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FIG. 9. Mass dependent analysis of the l" 1 + S and 
1 1 + P partial waves shows that the 2 waves are well de- 
scribed by two interfering Breit-Wigner resonances: the 1~~ 
state (a) fixed to the PDG parameters of the ^(1650) , 
which is the first observation of oj(1650) — > urj, and a new 
1+- state (b) with mass = 1594(±15) (+") MeV/c 2 and 
T = 384(±60)(+ 1 7 ° ) MeV/c 2 : c) shows the relative phase 
difference of l" 1 and 1 waves whereas d) exhibits the in- 
dividual Breit-Wigner phases and an overall constant produc- 
tion phase. 



^The calculation is mass dependent, and one expects the S 
to D wave ratio to increase for a lower mass h\ state due to 
angular momentum barrier effects. 



diets a J PC = 1+- hybrid near 1900-2000 MeV/c 2 @, 
and exotic states have been reported at masses lower than 
the flux-tube model expectations, an /ii(1595) would be 
a candidate for a hybrid or radial-hybrid mixed state. 



BACKGROUND STUDY 

The A distribution can be used to estimate the back- 
ground due to non-omega events under the omega sig- 
nal. The omega signal-to-background ratio estimated 
from Figure || is a little larger than 1 to 1. But since 
the omega is a spin 1 particle, the PWA uses the omega 
angular decay information to weight omega events more 
highly than non-omega events. 

To better understand the background under the omega 
and to see if these events are responsible for the ob- 
served structure in the l + ~0 + 5 and partial 
waves, events were selected from the sidebands around 
the omega region and a partial wave analysis of these 
events was performed. As expected, the \ S idebands dis- 
tribution, shown in Figure is flat and consistent with 
phase-space. 
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FIG. 10. Distribution of the sideband events as a function 
of the omega decay matrix element squared, A. Pure ui events 
would exhibit a distribution which linearly increases with A 
whereas phase-space events w/ould be flat. 



makes the background wave more favorable to describe 
those events which do not contain omegas. As seen in 
Figures Oa and [Tl|b, very few sideband events contribute 
to the F^O+S 1 and 1~~1 + P waves demonstrating that 
non-omega events are unlikely to be responsible for the 
observed resonant-like structures. The sideband events 
overwhelmingly contribute to the background intensity 
as shown in Figure 
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FIG. 11. The PWA results for the sidebands events. 
Shown are the : a) 1 1 + P intensity, b) V + S intensity, 
and c) the background intensity. 



A PWA of the sideband events was performed in the 
same manner as the original PWA of the data. Decay 
amplitudes were calculated using the sideband events 
whereas the normalization integrals and the list of partial 
waves included in the fit were identical to those used in 
the original PWA. Figure pi] shows results of the sideband 
PWA. Note that all partial waves except the background 
wave have an omega decay factor in the amplitude. This 



CONCLUSIONS 

We have collected a high statistics sample of the reac- 
tion 

at 18 GeV/c beam momentum, and performed a partial 
wave analysis on the cor] system where 



LO — > 7T 7T 7T 7T -^27 — ► 27. 

The lot] invariant mass distribution rises rapidly near 
threshold then increases slowly to a maximum near 1600 
MeV/c 2 . A PWA finds the data dominated by natural 
parity exchange waves. The results show 2 significant 
isoscalar J M e L partial wave intensities: a U + S 
intensity peaking near 1600 MeV/c 2 with a width of 
about 300 MeV/ c 2 , and a 1~~1 + _P intensity peaking near 
1650 MeV/c 2 also with a width of about 200 MeV/c 2 . 
In addition, a peak in the unnatural parity exchange 
wave 2 1~P exhibits a Breit-Wigner mass of 1830 ± 8 
MeV/c 2 and width of 86 ± 31MeV/c 2 . Unfortunately 
the phase of this state was not well determined and thus 
claim of a 2 state at this mass requires further con- 
firmation. Another notable feature was that the exotic 
partial waves for J PC = 0~~ and 2 H were found not to 
be required by the data. 

A mass dependent analysis of U + S and 1~~0 + P 
partial waves shows that the intensities and relative 
phase motion are well described by two resonating 
states. A 1~~1 + P state in good agreement with 
the PDG values for the w(1650), and a new /ii(1595) 
state with mass = 1594(±15)(+g°) MeV/c 2 and T = 
384(±60)(+ 1 7 ° ) MeV/c 2 . It is interesting to note that 
Close and Page predict that the isoscalar J PC = 1~~ 
hybrid should decay dominantly to states containing a 
vector meson such as pit and ior\ and not to traditional 
flux-tube decays such as bin. Possible interpretations for 
an /ii(1595) include an hi radial excitation, a hi hybrid, 
or a hi radial-hybrid mixture. 
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